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Emergent multicellular life cycles 
in filamentous bacteria owing 
to density-dependent 
population dynamics 
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Filamentous bacteria are the oldest and simplest known multicellular life forms. By using 
computer simulations and experiments that address cell division in a filamentous context, 
we investigate some of the ecological factors that can lead to the emergence of a multicellular 
life cycle in filamentous life forms. The model predicts that if cell division and death rates are 
dependent on the density of cells in a population, a predictable cycle between short and long 
filament lengths is produced. During exponential growth, there will be a predominance of 
multicellular filaments, while at carrying capacity, the population converges to a predomi- 
nance of short filaments and single cells. Model predictions are experimentally tested and 
confirmed in cultures of heterotrophic and phototrophic bacterial species. Furthermore, by 
developing a formulation of generation time in bacterial populations, it is shown that changes 
in generation time can alter length distributions. The theory predicts that given the same 
population growth curve and fitness, species with longer generation times have longer fila- 
ments during comparable population growth phases. Characterization of the environmental 
dependence of morphological properties such as length, and the number of cells per filament, 
helps in understanding the pre-existing conditions for the evolution of developmental cycles 
in simple multicellular organisms. Moreover, the theoretical prediction that strains with the 
same fitness can exhibit different lengths at comparable growth phases has important impli- 
cations. It demonstrates that differences in fitness attributed to morphology are not the sole 
explanation for the evolution of life cycles dominated by multicellularity. 
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1. INTRODUCTION 

Multicellularity is an organizational characteristic 
present in the majority of organisms whose size surpas- 
ses microbial scales. Phylogenetic inference suggests 
that evolutionary transitions to multicellularity have 
occurred several times during the history of life [1-6]. 
Nonetheless, despite its fundamental importance, it is 
difficult to empirically study the evolutionary and eco- 
logical forces that may lead to a transition from 
single-celled to multicellular organization. One set of 
theoretical explanations are based on examining the 
consequences of a shift between units of selection. 
These approaches consider the changing boundaries of 
an individual after a transition to multicellularity, and 
how the different activities of component cells can 
potentially lead to synergies that increase the fitness 
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of the multicellular assemblage [1,5,7,8]. Experimental 
results concerning social evolution in microbes are also 
largely compatible with this perspective [9-14]. 

An additional set of explanations for multicellula- 
rity are based on the potential selective advantages 
associated with increased size. Among the proposed 
size-related advantages are increases in the efficiencies 
of feeding [15-21], improved dispersal [16,22,23] and 
predator avoidance [24-28]. Although many of the pro- 
posed explanations are highly plausible, there is a lack 
of empirical validation. Some notable exceptions are 
measurements of motility and feeding efficiencies in 
Volvocalean green algae [20,22] and density-dependent 
growth in Myxoccoccus xanthus [17,21,29]. 

Multicellular bacteria with a filamentous form are 
ubiquitous in many environments, including freshwater, 
oceans, soil, extreme habitats and the human body 
[30-33]. Extensive empirical work has been done to 
examine and monitor filamentous bacteria that can be 
toxic or problematic in the environment [34-37]. Some 
filament-forming cyanobacteria also develop specialized 
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terminally differentiated cells, named heterocysts, that 
fix nitrogen and allow for the division of labour [38,39]. 
Although differentiation can represent a clear evolutio- 
nary advantage, theoretical and phylogenetic evidence 
suggests that in the cyanobacterial case, undifferentiated 
multicellularity evolved prior to differentiated multi- 
cellularity [40]. In aquatic bacteria, the most common 
hypothesis for the advantage of filamentation — and 
hence undifferentiated multicellularity — is an increase 
in size and the concomitant defence against predation 
by grazers [25-27,41-43] . However, experiments indicate 
that the avoidance of predatory grazers is not the only 
factor causing an increased frequency of filamentous bac- 
teria in aquatic environments [44]. Notably, in some 
bacterial species, filament formation appears to be depen- 
dent on the growth state of the population, whereby an 
increase in the dilution rate of chemostat cultures leads 
to longer filament lengths [27,45]. 

No theoretical studies address the distribution of fila- 
ment lengths and the population dynamics leading to 
shorter or longer filaments, although differences in 
length can reflect the extent to which a species is able 
to maintain multicellularity. Environmental conditions 
such as temperature, solar irradiation and nutrient 
concentrations have been detected as factors in deter- 
mining the mean size (filament length) of different 
cyanobacterial species [46-48]. Several of these factors 
also contribute to competition between species and adap- 
tation to different niches [49,50]. Filament breakage can 
occur because of external mechanical stress, lytic pro- 
cesses initiated by pathogens [51-53], or programmed 
cell death [54-58]. 

In the present study, we deliberately avoid a model- 
ling framework in which one assumes an a priori fitness 
advantage to multicellularity. We also do not aim to 
provide mechanistic explanations for the 'origin' of mul- 
ticellularity. Rather, by considering cellular growth and 
division in a filamentous context, we investigate the 
constraints that would affect filament formation and 
the maintenance of multicellularity. This gives an indi- 
cation of the capability of bacterial organisms to evolve 
multicellular life cycles according to their life-history 
traits. Our model is built on the basic assumption 
that the growth of a population of filamentous bacteria 
and the changes in the length of the filaments can be set 
in an ecological framework. In classical population 
dynamics, the change in population size is governed 
by the processes of birth and death (and sometimes 
migration), combined by the well-known logistic 
equation owing to Verhulst. In this model, birth and 
death rates are usually assumed to be decreasing and 
increasing functions of the population density, respect- 
ively. When the birth and death rates are the same, 
the population density is at a steady state and is said 
to be at its carrying capacity. The value at which the 
two rates are equal will be here referred to as the turn- 
over rate 6. The same carrying capacity can be achieved 
with different birth and death rate functions, and 
therefore at different corresponding turnover rates 
(figure la). 

Our approach to study distributions of bacterial fila- 
ment length places the described density-dependent 
concepts in the context of a population of cells that 



can form filamentous individuals. Although prokar- 
yotes had been commonly thought to be in principle 
immortal [59,60], there is now evidence for programmed 
cell death and lysis [61], and in addition ageing and 
senescence [62]. Taking this into account, a multicellu- 
lar filament is made up of several cells, and a population 
of filaments can also be considered as a population of 
cells. In this case, the population size is given by the 
total number of cells, while the birth and death rates 
coincide with the birth and death rates of each cell of 
the filaments. In such a population, one can vary the 
turnover at which a given fixed carrying capacity is 
achieved and ask the question if different turnovers cor- 
respond to different growth strategies; namely whether 
during its growth cycle the population is composed 
mostly of a few long filaments, or of many short fila- 
ments. Such a line of questioning is possible because 
as it will be shown, the rate at which a population 
reaches its carrying capacity can be different from the 
rate at which it reaches its stationary filament length 
distribution. The present work is an attempt to model 
and experimentally test the dynamics of filamentous bac- 
teria by only considering basic life-history traits of the 
population. The distribution of filament lengths in the 
different growth phases can provide insights into the 
developmental strategies of different filament-forming 
species, and the population conditions influencing 
multicellular development. 

In the following sections, we first state the model 
assumptions and we specify the experimental set-up 
for the culture of five filamentous bacterial species. 
Since turnover is a key factor in our model, we chose 
species that also exhibit different turnovers. We hence 
cultured heterotrophic species with fast generation 
time and high turnover, and photoautotrophic species 
characterized by a slower generation time and a low 
turnover. Heterotrophs differ from photoautotrophs in 
the way they obtain organic carbon for growth. While 
the former absorb organic compounds from the environ- 
ment, the latter autonomously fix carbon using solar 
light as an energy source. We present the simulation 
results and test the model through a comparison with 
our experimental results. 



2. METHODS 

2.1. Theoretical model 

We consider a population of multicellular, undifferen- 
tiated, non-branching filamentous bacteria, whose 
growth is regulated only by the density-dependent 
birth and death rates. In this context, the total 
number of cells of the population is used as the popu- 
lation size. The birth and death rates are represented 
by birth and death probabilities per iteration per-cell. 
At every iteration, each cell in a filament has a certain 
probability of dividing into two daughter cells or of 
dying. Based on these probabilities, the number of cell 
divisions and cell lyses are calculated in each filament. 
Accordingly, the filaments are elongated and then 
broken into smaller pieces. We simulate the change in 
time of the population until it has stabilized at its car- 
rying capacity. At each iteration, the following steps are 
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Figure 1. Simulation model (a) and results based on 1000 runs (b-d). (a) Schematic of birth and death rate functions dependent on total 
number of cells N c . The value of the population at which the two rates are the same is defined as the carrying capacity N c . We refer to the 
value of the rates at the carrying capacity as the turnover rate 8. The plot illustrates that the same carrying capacity can be achieved by 
species that have a different turnover (solid lines, species 1; dashed lines, species 2). (6) For all birth and death functions following the 
scheme in (a) N c reaches the steady state following the same growth curve despite the differing turnovers. The dashed line indicates 
the time at which the carrying capacity is reached (see mathematical definition of N c in §3.1.1). Error bars represent standard deviation 
of the mean, (c) When the population size reaches the carrying capacity (dashed line), the higher the turnover, the higher the total 
number of filaments Nf. Nonetheless, after many more iterations, N{ reaches a steady state of about 2500 filaments for all sufficiently 
small turnovers (plot not shown). Error bars represent standard deviation of number of filaments, (d) Box plots of filament lengths 
when N c reaches the carrying capacity. Populations with low turnover have higher median length values than populations with high turn- 
over. The corresponding pie charts on the top show the proportion of single-celled, double-celled and longer filaments in the population. 



performed (see illustration in the electronic supplemen- 
tary material, §S1). 

(i) Computation of the total number of cells in the 
population N c = X)j=i ^j'j where L is the maximal 
length of the filaments in the current generation 
and ^ is the number of filaments of length i. 

(ii) Computation of the birth and death proba- 
bilities per-cell per iteration, corresponding to 
the birth and death rate functions, respectively 
/3 = /3(JV c ) and 8= 8 (N c ): 



and 



S(N C ) = c 2 N c , c 2 = 



9 



(2.2) 



P(N 0 ) = - Cl N c +l, ci = 



1-0 



(2.1) 



where ^V c is the total number of cells and N c the 
carrying capacity. The latter is defined as 
the value of the population size, such that 
(3 (N* c ) = 8(N* C ) = 9 < 1, where 9 is the turnover. 
The birth rate is a decreasing function of N c . The 
death rate is an increasing function of iV c . 
(iii) For each filament, we compute the number of births 
b and number of deaths d occurring in the filament. 
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As /3 and 8 are the per-cell probabilities in a filament 
of length i, the total number of events b or d can be 
viewed as the number of successes in i trials, each one 
with probability of success f3 or 8. By picking from a 
binomial distribution B 1 we determine b and d for 
each filament of length i: 



and 



b = random pick from B(i, (3) 



d = random pick from B{ i, 8) . 



The filament is elongated by b cells. Then, a number 
of d cells computed on the basis of the original fila- 
ment length i are randomly selected for lysis in the 
elongated filament. The filaments resulting from 
the breakages are stored for the next generation. 
After completing the scan of all the filaments, the 
new population is set for the next iteration. The 
choice of initial conditions and the order at which 
birth and death are applied to filaments do not sig- 
nificantly affect the results of the simulations. The 
default initial condition for the simulations is a 
single cell. Alternative algorithms implementing a 
reverse or a random order of birth and death steps 
have been tested and produced similar results 
(figures not shown). 

As stated in point (ii), this model implements linear 
birth and death rate functions, with fixed carrying 
capacity. Alternatively, nonlinear functions have been 
implemented and provided qualitatively comparable 
results (see the electronic supplementary material, 
§S2). Considering the difference between those rates, 
given by f3 — 8 = — (c 1 + c 2 )N c + 1, one can observe 
that the net growth rate (c 1 + c 2 ) is independent of 
the turnover 8, as c 1 + c 2 = 1/N C . This implies that 
regardless of the turnover, the population size is 
expected to grow with the same curve for any turnover. 

2.2. Relation between turnover and generation 
time 

Before presenting the results, we provide a formula that 
helps to understand the connection between the model 
parameters and the experimental data. In the theoreti- 
cal model, the turnover is the characterizing property of 
a strain. However, in natural bacterial populations, 
turnover is not an easily measurable quantity. For bac- 
teria, doubling time during the exponential growth 
phase is usually calculated in lieu of generation time. 
Two populations with the same growth curve also 
have the same doubling time. Nonetheless, the reality 
is that the internal dynamics leading to this growth 
curve can vary. In the same time period, the same con- 
centration of bacteria can be achieved by undertaking 
many births and deaths, or by conversely having less 
births and also less deaths. This leads to different turn- 
overs and generation times. Moreover, at carrying 
capacity, the population density is stable, hence there 
is no indication of generation time at that stage by 
measuring doubling time. Here, we calculate the gener- 
ation time of bacterial populations at carrying capacity 
in relation to turnover rate. 



According to a common approach in ecology to 
describe the dynamics of an age-structured population 
[63], the mean age of parents at childbirth in a popu- 
lation with stable age distribution is calculated as: 



G = 



x— 0 ™ x 



(2.3) 



where x is the age, l x is the probability of survival from 
conception to age x and m x is the mean number of off- 
spring at age x. G can also be interpreted as a 
measure of generation time without making any 
assumption about the growth rate of the population 
[64]. In our case, age is measured in a number of iter- 
ations. We consider bacterial cells as individuals that 
give birth by cell division and die by a lytic process. 
Suppose that a cell survives a given iteration with 
probability 1 — 8. Then the probability that a cell is 
living after x iterations if 8 remains constant is given 
by lx~ (1 — S) x . If also the probability of cell division 
is f3 at every iteration, the average number of offspring 
at a given age x coincides with the probability of birth 
/3, whereby m x = /3. The generation time of a cell can 
then be expressed as: 



Godl(/3, 8) 



J2t 0 x(l-8y[3 = l-S 
E£=o(l-*r/8 8 ' 



(2.4) 



where the last equivalence was obtained by using geometric 
series. The standard formula in equation (2.3) applies when 
the growth rate is constant. Given the logistic nature of our 
model however, the functions (3 and 8 are actually depen- 
dent on the population size N c . Only when the population 
approaches the carrying capacity, namely when N c — ► N c , 
we have that j8, 8 — > 8, hence the growth rate is constant 
(equal to 0) . We then present a formula for the generation 
time at the carrying capacity. Substituting (5=8=8 
in equation (2.4), the generation time is given by 



G cc a(8, 8) 



(2.5) 



which is a decreasing function of the parameter 8 (graphi- 
cal illustrations in the electronic supplementary material, 
figure S2). Hence, a population with low turnover has a 
long generation time, whereby cell birth and deaths occur 
at a low rate. On the other hand, populations with high 
turnover have a short generation time, whereby births 
and deaths occur rapidly. 

2.3. Experimental setting 

2.3.1. Bacterial strains. The heterotrophic strains used 
in the experiments were Rudanella lutea DSM 19387 T 
[65] and two new bacteria isolated from a mud sample 
from tidal flats in Fedderwardersiel, on the North Sea 
coast of Germany, Fibrella aestuarina BUZ 2 T [66] and 
strain Fibrisoma limi BUZ 3 T [67]. The autotrophic 
strains were two axenic cyanobacteria, Nostoc muscorum 
SAG 25.82 (identical strains ATCC 27893; PCC 7120) 
and Anabaena variabilis SAG 1403-4b (identical strains 
CCAP 1403/4B; UTEX 377; ATCC 29211; PCC 6309). 

2.3.2. Batch culture experiments for heterotrophs. The 
growth of the three bacterial species was monitored in 
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batch cultures. Stock bacteria were grown on R2A 
plates for 3-4 days before starting the experiment. 
Colonies were homogenized in 0.7 per cent NaCl and 
finally distributed into 50 ml flasks containing 25 ml 
of SM (DSMZ 7) or R2A medium (DSMZ 830) for 
F. limi, R. lutea and F. aestuarina, respectively, reach- 
ing a starting optical density (OD) of 0.04. Bacteria 
were grown on a shaker (120g) at 29°C and after differ- 
ent time periods sub-samples were analysed for OD and 
filament length. For R. lutea and F. aestuarina., two 
batch cultures (A and B) were prepared and analysed, 
whereas for F. limi only one batch culture (A) was 
used for the experiment. Additionally, for each batch, 
three technical replicates (i.e. al, a2 and a3) were 
made. The OD was measured in a spectrophotometer 
(SpectraMax 384 Plus, Molecular Devices, USA) at 700 
and 600 nm. The samples for length measurement were 
fixed with formaldehyde (final concentration, 2%), 20 |jd 
were distributed on a microscope slide and pictures 
were taken with a digital camera (Color View, Soft 
Imaging System, 10X, 20X and 40X objective) con- 
nected to a phase-contrast microscope (Olympus BX 
51, Germany). For each time point, between 10 and 
20 pictures from separate fields were taken and the 
sizes of the filaments were individually determined by 
hand using the Soft Imaging software CellF (Olympus, 
Germany). At least 300 filaments from 10 different pic- 
tures were counted for each time point, in order to have 
a representative distribution of filament lengths. 

2.3.3. Batch culture experiments for cyanobacteria. As 
in the heterotrophic case, the growth of the two cyano- 
bacterial species was monitored in batch cultures. Stock 
cultures were grown in 20 ml of BG11 medium [68] in 50 
ml Erlenmeyers one week prior to the experiment. Two 
millilitres of grown culture were added to 18 ml of BG 11 
medium and cultivated on a shaker (150g) at room 
temperature and constant light (670 lux) for three 
weeks. Every 2 days, 200 (xl subsamples were taken, 
transferred into a microtitre plate and analysed for 
OD (at 430, 630 and 680 nm; SpectraMax 384 Plus). 
Filament length was measured by distributing 20 |xl 
samples in a counting chamber Neubauer-improved 
(Paul Marienfeld GmbH & Co, Germany) and pictures 
were taken with a digital camera connected to a phase- 
contrast microscope. Approximately 18 pictures with a 
10X magnification objective (or two pictures with a 4X 
magnification objective) were taken for each time point, 
and the sizes of the filament were determined using the 
Soft Imaging software CellF. For each time point, at 
least 50 filaments were counted in order to have a repre- 
sentative distribution of filament lengths. Three 
experiment replicates (Erlenmeyer flasks) of each cya- 
nobacteria and four measurement samples of each 
replicate were analysed for OD. Initial stock cultures 
were also examined. For filament lengths, four measure- 
ment samples were assayed from a single replicate. In 
figure 3d,e, we show the results of two experiment repli- 
cates (experiments A and B). 

To extrapolate the cell number from the filament 
length, an estimated size of cells forming a filament 
was measured from several light microscopy pictures. 



From at least 60 observations, the estimated mean 
cell lengths are 6.0, 6.7 and 5.9 |xm for R. lutea, 
F. aestuarina and F. limi, respectively. These values 
are higher than those measured from the inoculum of 
the same bacteria, indicating that cell size of hetero- 
trophs may change over time. For cyanobacteria, the 
mean cell length is less variable, and an approximate 
cell length of 4.0 |xm was used for the conversion. 

3. RESULTS 

3.1. Simulation results 

3.1.1. Turnover and fitness. For the rest of this article, 
we will refer to the time at which the population rea- 
ches its carrying capacity as the first time t, such that 
N c (t+ 1)/N c (t) <1 + e, where N c (t) is the population 
size at time t and e = 1%. The growth curve of the 
total number of cells gives information on the fitness 
of the strains. Figure lb shows that for any turnover, 
the total number of cells reaches the carrying capacity 
after less than 20 generations with the same growth 
curve (average and standard deviation of 1000 runs is 
shown). In our case, fitness is represented by net 
growth rate and does not depend on the turnover 
(§2.1). Hence, populations with different turnovers 
have the same fitness. On the other hand, when the 
total number of cells reaches the carrying capacity, 
the total number of filaments (iVf) differs significantly 
according to the turnover (figure lc). At this point, N{ 
ranges from values close to the carrying capacity, corre- 
sponding to almost unicellular filaments (8=0.9), to 
values under 20 units (0=0.001). Figure lb,c shows 
that the steady states of number of cells and number 
of filaments are not reached at the same time. This is 
because the stationary filament length distribution is 
reached much later than the steady state population 
density (carrying capacity). In the long run, the 
number of cells per filament reaches a stable equilibrium 
of about two cells per filament for sufficiently small 
turnovers, typically less than 0.1 (figure not shown). 
An analytical support for this result has been obtained 
by deriving an analogous model in continuous time, as 
presented in the electronic supplementary material, 
§S9. For higher turnovers, the effect of discretization 
leads to a deviation from this equilibrium. Hence, fila- 
ment length distributions can be affected by turnover 
rates in temporary phases of the population dynamics, 
before a stationary distribution is reached. An example 
is given by the time when the carrying capacity of the 
population is reached (17 iterations in this case). 

3.1.2. Length distribution at the time when carrying 
capacity is reached. Figure Id shows the box plot of fila- 
ment lengths when the total number of cells in the 
population (N c ) has reached the carrying capacity 
N c = 5000 cells. This time point however represents 
an intermediate stage before a stationary distribution 
of filament lengths is reached. The corresponding pie- 
charts indicate the fraction of single cells, double- 
celled filaments and filaments with more than two 
cells in the population. When turnover is close to 1 
(9= 0.9), the population consists mainly of unicellular 
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Figure 2. Sample simulations showing the mean filament length along successive iterations for two different turnovers, (a) 6 = 
0.9; (6) 6 = 0.1. Squares indicate the mean of the filament lengths obtained through 1000 runs. Edges of the error bars correspond 
to the 2.5th and 97.5th percentiles of the distribution of the mean filament lengths in the 1000 runs. In both cases, the mean 
length reaches a peak in the transient phase, and then decreases gradually towards a state comparable to the one of the beginning. 
Lower turnovers allow for higher transient peaks. The dotted line indicates the iteration at which the carrying capacity is reached. 



and double-celled filaments. Note that 0 = 0.9 means 
that at carrying capacity, each cell has a 90 per cent 
chance of dividing as well as a 90 per cent chance of sub- 
sequently dying. When 0 = 0.1, filaments with one or 
two cells represent about half of the population. Fol- 
lower turnovers (0=0.01, 0.001), the median values 
of the length are significantly higher. In the latter 
cases, filaments with more than two cells are more 
than 90 per cent. Figure Id shows that, as the turnover 
increases, the fraction of short filaments when the popu- 
lation reaches its carrying capacity also increases. The 
distribution means at this growth stage were statisti- 
cally compared by means of a one-way ANOVA and a 
multiple comparison (see the electronic supplementary 
material, §S6). The means obtained for 0=0.9 and 
0=0.1 were not significantly different. All the other 
pairwise comparisons were significant. 

3.1.3. Length distribution during the transient phase. 
Figure 2 shows two examples where the mean length 
of the filaments was calculated every iteration. The 
plot shows the average of 1000 runs, with error bars 
indicating the 2.5th and 97.5th percentiles for the 
mean. The dotted line indicates the iteration at which 
the total number of cells reaches the carrying capacity. 
An increase to a peak of the mean length followed by a 
decay towards a steady state is observed in all cases. 
At the end of the decreasing phase, the mean length 
is comparable to that in the beginning. The same pat- 
tern was also observed for 0=0.01, 0.001 (see the 
electronic supplementary material, §S3). An analytical 
characterization of the filament length dynamics has 
been carried out in the electronic supplementary 
material (§S8), where an approximation for the iteration 
of first filament breakage has been derived. In order to 
statistically analyse the average length trend, we con- 
sidered three reference time points and we compared 
the distribution of the means at those points using a 
non-parametric hypothesis test. As reference points we 
chose the starting point, the peak and the carrying 
capacity. The outcome of the tests provides statistical 



evidence of a significant increase and subsequent 
decrease of the average length. Details and results are 
given in the electronic supplementary material, §S6. 

Significant differences are found in the maximum 
average length achieved. Long filaments, with a mean 
length up to 250 cells, are observed at a turnover of 0 = 
0.1. Longer filaments can be achieved with lower turn- 
overs (electronic supplementary material, figure S4). 
For a higher turnover ( 0 = 0.9) , the average length of fila- 
ments at its peak does not exceed 100 cells. The observed 
trend indicates the prevalence of a cycle in the mean 
length induced by the birth and death rates, whereby 
filaments are short in the beginning, longer in the transi- 
ent phase and again short at the carrying capacity. 
Moreover, in the transient phase, filaments with low 
turnover can elongate more than filaments with high 
turnover. 

3.2. Experimental results and model validation 

In the theoretical model, the population stays at the 
carrying capacity once this has been reached. However, 
as bacterial strains were not cultured in a chemostat, 
carrying capacities could not be maintained indefi- 
nitely. In order to assure that the experimental 
dataset and the simulation data are comparable, we 
indicate with a dotted line the time point where the 
bacterial populations reach their carrying capacity on 
each plot of the experimental results. This point was 
determined from the OD curves (not shown). Sub- 
sequent to reaching the carrying capacity, the process 
of filament breakage is eventually accelerated in the 
empirical case owing to a progressive decrease in the 
available nutrients. 

In order to prove the significant increase and 
decrease of bacterial filament length, we performed an 
ANOVA test associated with a multiple comparison 
on the data based on Tukey's HSD test [69]. Details 
of the tests for each species and within species are pro- 
vided in the electronic supplementary material, §S6. 

In the following sections, the time point 0 is rep- 
resented by the inoculum. The starting culture of 
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Figure 3. (a) -(e) Empirical measurements showing mean length of bacteria against time. Edges of the error bars correspond to 
the 2.5th and 97.5th percentiles of the measured lengths. Diamonds indicate the mean of measured filament length at each time 
point. Data from both heterotrophic (a- c) and photoautotrophic (d- e) species are shown. The zero time point corresponds to the 
inoculum. The dotted line indicates the time at which the carrying capacity was reached. (/) Box plots of filament lengths when 
the bacterial populations reach their carrying capacity. The capital letter next to the species name indicates the experiment used 
in the plot. The corresponding pie charts on the top indicate the proportion of single-celled (white), two-celled (grey) and longer 
(black) filaments in the population. 



heterotrophs was taken from plate colonies and vor- 
texed. Because of the nutrient depletion in plate 
colonies and mechanical breakage through vortexing, 
the heterotrophic inoculum is composed mainly of 
single cells. The starting cultures of cyanobacteria 
were taken from cultures grown for more than one 
week, gently pipetted and diluted. At inoculum, 



cyanobacteria filaments had an average filament 
length of about 45 cells (N. muscorum) and about 75 
cells (A. variabilis). 

3.2.1. Mean length trend. Figure 3 a- c show the change 
in time of the mean length of heterotrophic bacteria 
R. lutea, F. limi and F. aestuarina. The pattern is 
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Figure 4. Bright field micrographs of F. aestuarina bacteria at different stages of their growth, (a) The initial population 
(inoculum) is composed of many small filaments, (b) After 17 h, the filaments reach their maximal length, (c) After this 
peak, the bacteria progressively break towards a mixed population of short and medium filaments as they reach their carry- 
ing capacity (48 h). (d) Eventually, the filament lengths approach the mean of the inoculum (168 h). Scale bars: (a,d) 50 |xm; 
(6,c) lOOujn. 



similar in all three cases: the mean length increases until 
a peak, after which it decreases and reaches a value 
approaching that of the inoculum (see the electronic 
supplementary material, table SI). This trend is 
recognizable in figure 4, showing micrographs of 
F. aestuarina taken at different time points. Remark- 
ably, experiments A and B on F. aestuarina give very 
similar results, showing a high repeatability of the 
observed pattern (figure 3c). A comparison of the 
theoretical and F. aestuarina length distributions 
against time is shown in the electronic supplementary 
material (§S4 and figure S6). The mean length of cyano- 
bacteria, tracked for two months, shows a qualitative 
trend similar to that of heterotrophs, albeit much 
slower. Both the A. variabilis and N. muscorum species 
have a mean-length peak in the transient phase. How- 
ever, because of the low turnover, they need more 
time than heterotrophs to reach the stationary distri- 
bution, as predicted by the model. When they reach 
the carrying capacity, they are still long. 

The pattern observed in heterotrophs is comparable 
with the simulation results shown in figure 2. The 
trend of the cyanobacteria, although not so markedly 
conclusive, gives an indication that the model pre- 
dictions also hold in autotrophs. An additional 
illustration of the mean length trend of the cultured 



species is provided in the electronic supplementary 
material, figure S5. 



3.2.2. Length in the transient phase. Generation times in 
heterotrophs and cyanobacteria are very different. 
While the former are observed on the scale of hours, 
the latter grow on a scale of days/weeks. The difference 
could be, in principle, owing to the different trophic 
state of the bacteria. This hypothesis could be tested 
by means of a comparison between different generation 
times intraspecifically in the autotrophs or in the het- 
erotrophs. However, this goes beyond the scope of the 
present work and we proceed relying on the mathemat- 
ical derivation of §2.2. According to the relationship 
between generation time and turnover established in 
§2.2, we assume that the heterotrophs and cyanobac- 
teria have high and low turnovers, respectively. The 
simulation results presented in §3.1.3 show that the 
peak value of the mean length in the transient phase 
increases with decreasing turnover. Figure 3 shows 
that heterotrophic filaments will rarely be longer than 
60 cells. The length of cyanobacteria can be significantly 
higher than that of heterotrophs. Nostoc muscorum and 
A. variabilis filaments have up to around 350 and 600 
cells, respectively. Hence, cultured bacterial populations 
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the trend of each culture, (a) Computer simulations of successive transfers. When the average filament length is at its peak, the 
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support the model predictions, indicating that bacteria 
with a putatively higher turnover (in this case hetero- 
trophic species) are able to elongate less than those with 
lower turnover (cyanobacteria). 

3.2.3. Distribution of lengths at the time when carrying 
capacity is reached. In figure 3/, the box plots and the 
pie charts show the proportion of short and long fila- 
ments in the cultured bacterial populations when they 
reach their carrying capacity. As for the simulation 
results, the pie charts illustrate three length classes, 
namely single cells, double-celled filaments and 
filaments with more than two cells. The heterotrophs 
show a considerable amount of short filaments, with a pro- 
portion of single cells and double-celled filaments varying 
across species in the range 8-16% and 29-42%, respect- 
ively. In cyanobacteria, filaments with one or two cells 
are absent. Figure 3/ can be compared with figure Id. 
Again, the experimental results confirm the theoretical 
predictions, according to which short filaments at the car- 
rying capacity are more abundant in populations with 
high turnover. 



3.2.4- Serial transfer experiment and simulations. If fila- 
ment elongation is governed by density-dependent 
processes as hypothesized in the model in §2.1, then 
interventions on such underlying processes should pro- 
duce testable hypotheses. Experimental manipulation 
of cell densities is one simple possibility. If the popu- 
lation drops significantly below its carrying capacity, 
the filaments that survive in this new condition 
should be able to elongate or at least maintain their 
original length until the carrying capacity is reached 
again. We simulated and experimentally carried out a 
series of successive transfers of bacterial populations, 
whereby a fraction of the population was transferred 
to a fresh medium whenever the average length of the 
original population was close to its peak. Filaments 
developing in the fresh medium were expected to keep 
or increase the length reached at the previous stage. 
Details of the computational and experimental settings 
are provided in the electronic supplementary material, 
§S7. Figure 5a shows the average length of filaments 
obtained by successive transfer simulations. Figure 56 
shows the average length of F. aestuarina tracked 
along successive transfers. In this case, filament size 
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slightly increases or remains almost constant in most of 
the cases. The fact that filaments do not elongate so 
strongly as in the simulations could be owing to physical 
handling (breakage of long filaments while pipetting) or 
to other culture conditions influencing the elongation 
that were not considered in the model. The model 
results indicate that the filaments transferred to the 
fresh medium maintain or increase their length, and 
eventually, start shortening when they reach their car- 
rying capacity again. The experimental results, 
although not as strong as the simulations, show that 
the transferred bacteria are generally able to maintain 
their length. 



4. DISCUSSION 

Filamentous bacteria can be viewed as a simple example 
of multicellular individuals. Some of them, such as fila- 
mentous cyanobacteria, have been present in the 
environment since the early evolution of life, as 
indicated by fossil records and phylogenetic analyses 
[70-73]. The theoretical model shows that at various 
life stages, two populations with the same fitness and 
carrying capacity can differ in the composition of fila- 
ment length classes according to the life-history traits 
regulating their growth. Although the long-term 
stationary distribution of filament length is the same 
for any sufficiently small turnover rate, the differences 
owing to turnover observed at other growth stages 
(e.g. transient phase, when reaching carrying capacity) 
are quantitatively significant and biologically relevant. 
Filaments belonging to a population with high turnover 
(and hence short generation time) have moderate 
lengths at their maximum growth phase. Furthermore, 
when the population reaches carrying capacity, they are 
chiefly unicellular. Individuals of a population with low 
turnover (and hence long generation time) can instead 
become very long in the transient phase, and still be 
multicellular when reaching the carrying capacity. The 
species that we cultured show evidence that the predic- 
tions of the model hold. The relationship between 
bacterial generation time and turnover derived in §2.2 
helps us to understand the relationship between simu- 
lations and experiments. The heterotrophs have a faster 
generation time than the cyanobacteria (and hence a 
higher turnover). As a consequence, while the hetero- 
trophic filaments are rarely longer than 60 cells, the 
cyanobacteria can reach lengths of up to 600 cells. Both 
theory and experiments show that according to different 
life-history traits, bacteria can cover a wide spectrum of 
lengths, indicating that multicellularity can be achieved 
to different degrees. Moreover, a key feature of this 
model is that the predominance of multicellularity 
during the life cycle is not necessarily determined by mor- 
phology-dependent differences in fitness. Here, the factors 
affecting multicellularity are birth and death rates, and 
hence turnover rate and generation time. 

The model indicates a common temporal behaviour 
for filament length distributions. The mean length 
increases and peaks during the transient growth 
phase, and then decreases until it reaches a value 
approaching the initial conditions. The results from 



cultured bacterial populations of heterotrophic and 
photoautotrophic species show that the predicted 
trend in filament length distribution is robust across 
species, and is representative of what occurs in nature. 
The fact that individuals eventually become shorter 
(and of average length 2 at the stationary distribution) 
regardless of the turnover is because of the effect of cell 
lysis on a filament. Consider a filament of length L. 
A cell division (birth) increases its length by one unit, 
namely the filament length becomes L+l. By contrast, 
after a cell lysis, the expected filament length will be on 
average L/2. The average effect of cell death on length 
is hence much larger than that of cell birth. This is con- 
sistent with the expectation that the average filament 
length in a population can keep increasing if the death 
rate remains sufficiently low, as shown in the estimation 
of the time of first filament split (electronic supplemen- 
tary material, §8). When the average filament length 
in the population is at its stationary distribution, cell 
birth and death have the same effect. This is achievable 
only if the filaments are of average length 2. This argu- 
ment is supported by the analysis of the analogous 
continuous time model derived in the electronic sup- 
plementary material. In this context, a way to avoid 
the breakage into small filaments is to decrease the 
effect of cell lysis on organism size. A solution in this 
direction is the evolution of a second or third dimension, 
as in spherical multicellular bacteria [74] or microbial 
biofilms. In this case, if we think of a square or a 
ball of N cells, a cell division will still lead to an 
individual oi N+ 1 cells, but cell death will also only 
decrease size by one unit, resulting in a size of N — I 
(as opposed to N/2). 

There are, to our knowledge, no previous theoretical 
and experimental results that are directly comparable 
to the results presented here. However, the extensive 
work of Hahn & Hoefle [27] and Hahn et al. [45] set 
important precedents and reference points. These 
works focus on the effects of predators on selection for 
longer filament lengths. Their results are not directly 
comparable to those presented here because the death 
rates owing to predators preferentially afflict smaller 
filaments and single cells. However, they also provide 
revealing information on the dynamics of filament for- 
mation in the absence of predators. Hahn et al. [45] 
show that increasing the nutrient flow and the dilution 
rate of a chemostat in the absence of predators leads to 
longer filament lengths. This corresponds to our predic- 
tion that in transient growth conditions with lower 
population densities, filament lengths would get 
longer. Furthermore, although they do not present 
extensive data on this, [45] they mention that in 
batch culture, the percentage of multicellular filaments 
is higher in the exponential phase in comparison with 
the stationary phase. This situation is again consistent 
with the predictions presented here. 

The model studied here is a general model of filamen- 
tous growth, and is not restricted to prokaryotes. Many 
simple forms of eukaryotic algae grow in a multicellular 
filamentous form, and the theoretical results presented 
here are in principle also valid for such eukaryotes. 
However, it is clear that eukaryotic algae have sur- 
passed their prokaryotic counterparts in their 
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morphological diversity. There are at least six major 
algal lineages in which filamentous multicellularity 
has evolved: the Phaeophyta, Rhodophyta, Chryso- 
phyta, Chlorophyta, Charophyta and Embryophyta 
[75]. However, in each of the latter cases, higher dimen- 
sional variants have also evolved, whereupon cell 
division can occur in more than one planar dimension. 
This diversity can for example range from the simple 
globular multicellularity of some of the Volvocales 
(within the Cholorphyta), to the highly differentiated 
and developmentally complex vascular land plants 
(within the Embryophyta). Furthermore, plants can 
differentiate into diverse developmental modules such 
as branching nodes, leaves, inflorescences and root 
structures. By combining such modules in different 
combinations, they can access a wider variety of mor- 
phologies during their growth process [76,77]. This 
flexibility in morphology can translate into both 
phenotypic plasticity and evolvability. 

Developmental life cycles are not unique to multicel- 
lular eukaryotes, and an understanding of how they 
evolved is lacking. Theory and experiments in this 
article indicate that the size of the individual, in 
terms of number of cells, can follow a cycle: from unicel- 
lular (or short length) to multicellular and back to 
unicellular (or short length) again. Each time, the 
initiation of the cycle can be hypothetically achieved 
by moving cell density below the carrying capacity 
(e.g. increases in fluid or nutrient availability). The 
transfer experiment and simulations do not contradict 
this hypothesis. 

In our model, only the birth and death rates play a 
role in the change in size during the bacterial life 
cycle. These life-history traits are intrinsic properties 
of every living organism. This pattern, which automati- 
cally arises from the interplay between ecology and the 
filamentous nature of the bacteria, is an emergent prop- 
erty. As such, it can serve as the basis for a primitive life 
cycle, upon which a more complex developmental pro- 
gramme can be subsequently built. As a case in point, 
some species of Nostoc can differentiate into hormo- 
gonia during low growth phases characterized by high 
stress or scarce resources. These hormogonia are short 
filaments specialized for survival in harsh conditions, 
which subsequently migrate to a new area before growing 
into full-grown filaments again [78] . Our results suggest 
that the cycling of filament lengths between growth 
and static phases of the population is a pre-existing con- 
text within which multicellular developmental cycles 
and differentiation can evolve. 
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